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Abstract 

We review equilibrium thermodynamic properties of systems of magnetic particles like ferrofluids 
in which dipolar interactions play an important role. The review is focussed on two subjects: (i) the 
magnetization with the initial magnetic susceptibility as a special case and (ii) the phase transition 
behavior. Here the condensation ("gas/liquid") transition in the subsystem of the suspended 
particles is treated as well as the isotropic/ferromagnetic transition to a state with spontaneously 
generated long-range magnetic order. 
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I. INTRODUCTION 

Ferrofluids Q| are suspensions of ferromagnetic or ferrimagnetic particles with diameters 
of the order of 10 nm in a liquid carrier. A common combination is magnetite (FeaC^) in 
hydrocarbons. The volume fraction of the magnetic material is typically a few percent. The 
particles are monodomain. They are too small for the build-up of domain walls . Thus, 
the quantum-mechanical exchange interaction between the electronic spins can keep all the 
atomic magnetic moments aligned so that the total magnetic moment, m, of the particles is 
of the order of about 10 4 ^b- Hence, to describe the equilibrium magnetic properties of the 
ferroffuid one can take each particle i to have a magnetic moment nij(t) that fluctuates in 
time, however, with its magnitude m; = |rrij(t)| being fixed. Thus, we do not address here the 
magnetic properties of the ferrofluid system on the nm scale, i.e., of the magnetic particles 
themselves: Questions related to the size and form of the magnetocrystalline anisotropy, to 
magnetostrictive effects, to shape anisotropy, to dead layers etc. of the particles n □ □ □ 
are not addressed here. 

Whether the magnetic moment mj(t) fluctuates only as a result of the rotational diffusion 
of the particle in the carrier liquid (Brownian relaxation) or also because the temperature 
is sufficiently high to stimulate rnj(t) itself to move in the potential landscape over the 
anisotropy energy barriers of the particle's lattice structure (Neel relaxation) does not influ- 
ence the long-time equilibrium magnetic properties of the suspension but only its dynamics. 
In any case the magnetic particles cause the ferrofluid to display the equilibrium thermody- 
namics of a paramagnetic liquid — if one postpones for the moment the currently debated 
question of whether spontaneous long-range magnetic order is possible or not. The initial 
susceptibility can easily reach values of x ~ 1 an d higher. Being orders of magnitude larger 
than for ordinary (molecular) paramagnetic liquids one thus speaks of a "superparamag- 
netic" liquid. 

Ferrofluid particles are small enough to avoid segregation caused by gravity or customary 
magnetic gradient fields. But the strong van der Waals attraction between the nanoscale 
particles would lead to irreversible agglomeration. To prevent this, they are either coated 
with polymer surfactants of about 2 nm length whose entropic repulsion counteracts the van 
der Waals forces, or with ionic groups that prevent agglomeration due to their electrostatic 
repulsion. In the latter case the carrier is a polar liquid like water, containing the necessary 
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counterions. 

In a sufficiently diluted ferrofluid the magnetic particles can be thought as noninteracting, 
and the equilibrium magnetic properties of such a ferrofluid are those of an ideal param- 
agnetic gas. Still, there is an important difference to a molecular paramagnetic fluid like, 
say, oxygen: the ferrofluid particles are not identical, they differ both in size and magnetic 
moment. This polydispersity has influence on the properties of the ferrofluid, in particular 
also on the equilibrium magnetization. 

At higher concentrations the ideal gas approximation fails to reproduce the equilibrium 
magnetization curve of ferrofluids. This is already the case for susceptibilities of the order 
of unity. Then, the particle interaction — and first of all the dipolar interaction — has to 
be taken into account. Since there exists a perfect analogy between magnetic and electric 
dipoles, the large number of theories developed for the latter case are directly applicable to 
this problem. 

The thermodynamics of dipolar interacting particles is not only of interest because of 
the applicability to (electric) polar fluids and ferrofluids. It is also of large theoretical 
interest because the dipole-dipole interaction shows at least two interesting features: (i) it 
is attractive and/or repulsive, depending on the orientation of the particles and (ii) it is of 
long-range nature. 

Because the dipolar interaction is neither always attractive nor always repulsive it is not 
trivially clear how the magnetic moments of the particles affect, e. g., the condensation phase 
transition of the suspended particles and especially whether the dipole-dipole interaction 
can cause such a transition in a system of particles without additional attractive interactions. 
This question is still being discussed in the literature. An externally applied magnetic field 
can be expected to further modify the phase transition behavior. 

Another example of a phase transition where the magnetic moments play an important 
role is the isotropic/ferromagnetic transition, i.e., a spontaneous breaking of the rotational 
symmetry resulting in a nonvanishing magnetization without applied magnetic field. Also 
this question whether the dipolar interaction alone can trigger such a phase transition is not 
yet settled. 

There is quite an extensive literature on ferrofluids. The references to it are accumulated 
in the special issues of the Journal of Magnetism and Magnetic Materials that cover the 
International Conferences on Magnetic Fluids. Odenbach et al. p show in the latest 
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issue that by 2002 the total number of papers on ferrofluids has grown to well above 6000 
with about 1000 being published in the period 1998-2001 [7J. The research topics cover 
problems as different as fluid dynamics and the modifications in classical pattern formation 
experiments by using ferrofluids, pattern formation processes genuine to ferrofluids as the 
Rosensweig instability or the labyrinth instability, magnetorheological prop erties, dynamical 



We restrict ourselves here to the topics where the dipolar interaction between the magnetic 
particles is of major importance, namely the already mentioned equilibrium magnetization 
and initial susceptibility on the one hand, and phase transitions on the other hand. 

The review is organized as follows. Following this introduction we discuss in Sec. HT1 basic 
statistical mechanics of systems of dipolar particles. In Sec. IHII the influence of the dipolar 
interaction on the initial susceptibility of systems of otherwise spherical symmetric particles 
is reviewed, mainly for dipolar hard spheres and Lennard- Jones particles with an addi- 
tional magnetic moment (Stockmayer particles). This question has already been discussed 
for the electrical case before ferrofluids became an active research area. The more general 
question of the equilibrium magnetization on the other hand was mainly investigated with 
regard to the application to ferrofluids. The literature concerning the equilibrium magne- 
tization is reviewed in Sec. IIVI In Sec. |V] the effect of additional nonspherical interactions 
other than the dipolar one are briefly discussed. Phase transitions of the gas/liquid and 
isotropic/ferro magnetic type are reviewed in Sec. I VII We conclude in Sec. IVHI 

Because of the analogy between systems of magnetically and electrically polar particles 
and for reasons of simplicity, we will always use the magnetical terminology when reviewing 
the literature, even when only the electric case is considered. This is mostly the case in the 
older literature. 



magnetic susceptibility, medical and engineering applications 
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II. STATISTICAL MECHANICAL DESCRIPTION 

A complete description of ferrofluids within the framework of statistical mechanics would 
have to include all the components of a ferrofluid: magnetic particles, carrier liquid, polymer 
surfactants etc. But theoretical approaches to explain the thermodynamic properties con- 
nected to magnetism generally concentrate quite successfully on the magnetic component. 
Since the magnetic properties of the carrier liquid can be neglected, it is clear that it can 
only have indirect influence on the magnetization. A more elaborate model of a ferrofluid 
describing both the magnetic particles and the carrier fluid is discussed by Kalikmanov jl^ , 
who shows that under some assumptions the carrier fluid has no influence on the equilibrium 
magnetic properties of the ferrofluid whole. 

Ignoring the carrier, the ferrofluid can be described as a system of particles, i = 1 ...N, 
each carrying a magnetic moment rrij and thus interacting with an external magnetic field 
H e via the potential 

Vi = -mi ■ H e . (2.1) 

The magnetic moments as used here include the vacuum permeability /io, see e.g. Eq. (12. 4j) . 
The absolute values |rrij| are fixed; but in general they are different from each other. These 
magnetic particles interact pairwise with each other via the dipole-dipole interaction poten- 
tial 

Here = r^r^ is the distance vector between particles i and j (see Fig. |TJ). Since this 
potential decreases radially only as r~ 3 the dipolar interaction is of long-range nature. 

In addition to ()2.2j) there generally exist also some short-range interactions between the 
particles that we capture by the potential VM R . It might include nonmagnetic components 
and in general also short-range magnetic interactions involving quadrupolar or higher order 
magnetic moments. For later use we define the dimensionless potentials 

y. V DD V SR 

y% DD 13 SR «? rrt o\ 

The long-range nature of the magnetic interaction requires some care in statistical me- 
chanical considerations. The local magnetic field acting on a particle is the sum of the 
external field and the dipolar fields of the other particles. The dipolar contribution depends 
on the distribution of the particles on a macroscopic scale, i. e. in particular also on the 
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probe geometry. In macroscopic magnetostatics this property translates into the fact that 
the equilibrium magnetization M of a magnetic medium depends on the internal field H. 
This macroscopic field depends itself on H e , M, and geometry. The difference between H 
and H e can therefore be interpreted as an effect of dipolar interaction. 

A. Ideal paramagnetism 

If all dipolar interactions are discarded then the magnetic particles of the ferrofluid feel 
only the external magnetic field, i.e., H = H e . This is an useful approximation if M -C H e 
such that H rs H e . Assuming in addition all particles to be identical - the monodisperse 
case - with common magnetic moment |m.^| = m, then the equilibrium magnetization is 
given by 



with H = H e . Here N/V is the particle density and £ is the Langevin function, C(a) = 
coth(a) — 1/a. Its argument a = mH/ k-^T measures the energy of the moment m in the field 
H = H e in units of the thermal energy k^T . The saturation magnetization M sat = mN/ 
of the ferrofluid corresponds to the case of parallel dipoles. 

One might think that using in Eq. ()2.4)1 the real internal field H in the ferrofluid instead 
of H e would be an improvement of the ideal paramagnetic model by itself. However such 
a replacement does not incorporate the effect of dipolar interactions in a systematic and 
controlled manner since they additionally modify the simple functional relation ()2.4|) between 
equilibrium magnetization M and internal field H in various ways that are of the same order 
as the replacement of H e by H in Eq. Ij'2.4j) . 

In a polydisperse extension of the above described ideal, noninteracting monodisperse 
model the magnetic moments of the particles are allowed to be different. Then the magne- 
tization is given by 



Here the sum extends over all magnetic particles, % = 1,...,N. Using the mean magnetic 
moment m = (1/A) m % to introduce the reduced moment /ij = rrii/rn of particle i we can 
write more conveniently 




(2.4) 




(2.5) 




(2.6) 
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where 

M sat = ^ , £ poly (5) = ^ E ^ £ ^ ' 5 = j^f ■ ( 2 ' ? ) 
B. Weiss and Onsager model 

The earliest model of a selfinteracting magnetic medium is the mean-field Weiss model 
Q|. In this monodisperse model every particle is thought to be located in the center of an 
empty spherical cavity that is surrounded by a magnetic continuum with an internal field 
H and magnetization M. In thermal equilibrium, the magnetization is given by a Langevin 
function into which enters the local field Hi oca i = H + M/3 within the spherical cavity. This 
leads to the implicit relation 

m ( M\ 
H + —\ 



m = m -> c [w{ h+ t)\ • < 2 - 8) 

that can be solved numerically to give the sought after equilibrium magnetization M(H). 
The initial susceptibility of the Weiss model 

X = ^(H = 0)= XL - , (2.9) 

is equivalent to the well known Clausius-Mosotti relation. Here 

K-agv-v-^ ( 2 ,o) 

is the Langevin initial susceptibility of the ideal monodisperse paramagnetism ([2.4)1 . 

The Weiss model works well for weakly interacting ferrofluids but strongly overestimates 
the magnetization of concentrated magnetic fluids. It predicts ferromagnetic solutions, i. e., 
a spontaneous magnetization M(H = 0) ^ whenever xl > 3. While this criterion is at 
odds with experiments the possible existence of ferromagnetic phases in systems of dipolar 
hard spheres cannot be ruled out for stronger interactions (c.f. Sec. IVI B|) . 

A similar early approach to the problem of a selfinteracting magnetic medium is the 
Onsager theory [l5| originally conceived for polarizable molecules. Therein, say, spherical 
molecules occupy a cavity in polarizable continuum with susceptibility x- The field acting 
on the molecule is the sum of a cavity field plus a reaction field that is parallel to the actual 
total (permanent and induced) moment of the molecule. The resulting equation between xl 
and x reads here 

X(3 + 2 X ) ,„„. 

WTxf = XL - (2 - n) 
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While the Weiss model overestimates the initial susceptibility of concentrated ferrofluids, 
the Onsager model generally underestimates it. It is instructive to write \ according to (|2.9J1 
and (j2.11j) for small \l as power series in \l 

X = Xl + \xI±\xI + - ■ (2.12) 

Here the positive sign refers to the Weiss model and the negative one to the Onsager model. 
The first and second term agree with a systematic expansion of x in terms of dipolar coupling 
strength and magnetic volume fraction (c.f. Sec. IIII A|) . However, the third order term of 
this expansion reads +^Xl ( see Eq. 13.8)1 . Thus both, Onsager as well as Weiss model 
strongly overestimate the influence of this term. 



C. Short range potentials 



If one wants to go beyond the mean field models of Sec. Ill Bl one has to specify the 
potential with which the magnetic particles interact besides the long-range dipolar one. 

Rosensweig [1] has proposed a quite elaborate potential for the short-range interaction, 
V SR , of sterically stabilized magnetic particles. It consists of the sum of the van der Waals 
attraction between spherical magnetic particles and the sterical repulsion of their surfactants, 
i.e., V = V vdW + V ster . The van der Waals attraction between the spheres is given by 



VZ 



vdW, 



A 
"6 



+ 



P + Al (Z + 2) 2 
Here A is the Hamaker constant. The quantity 

2 S 2 7™a q 



+ In 



I 2 + 41 
{l + 2f 



(2.13) 



D r 



D r 



(2.14) 



is given by twice the surface-to-surface distance s of the particles reduced by the diameter 
-D core of the magnetic particles while is the distance between their centers. The sterical 
repulsion of the surfactants is described by 



■^ster , 



TcD*£kBT 



2- l A 
t 



l + t 
1 + //2 



< / < 2t . 



(2.15) 



Here £ is the surface density of the polymers and 

2<5polymer 



t 



D r 



[2.16) 
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is given by twice the length, S po i ymer , of the surfactant polymers. In the Rosensweig potential 
the carrier liquid component is not completely ignored but rather taken into account on a 
continuum level via the electric permeability of the liquid that enters into the Hamaker 
constant A. 

Fig. shows the interaction potential V SR = V vdW + V stcr (pTHH E35D for typical values 
given by Rosensweig: A = 10~ 19 J, £ = lnm~ 2 , D core = 10 nm, and <5 po iymer = 2nm. The 
need for sterical stabilization becomes apparent in this plot. The van der Waals interaction 
is strongly attractive for small particle separations, it even diverges in this simple model. 
The sterical potential modifies the interaction for distances less than 4 nm, such that only 
a small attractive tail at larger distances remains in the combined potential. At very small 
surface-to-surface separations the potential still diverges, but this inner attractive region is 
effectively shielded by a potential barrier with a height of about 25 k^T. 

To obtain an equivalent model potential for electrostatically stabilized ferrofluids, V ster 
is replaced (see e. g. by the interaction potential 



of charged spheres in ionic solutions. Here is the Debye length, e is the dielectric constant 
of the carrier liquid, and Q e fi is an effective charge of the particles. The combined potential 
shows qualitatively the same features as in the sterical case. 

Another important model potential for describing dipolar particles is the Stockmayer 
potential in which the short-range interaction is given by the Lennard- Jones function for 
the van der Waals interaction of point-like particles 



Here D is the collision parameter marked by the zero of this potential and Vq is the depth 
of the potential well. 

However, most often the effects of particle interaction are discussed within the simpler 
model of dipolar hard spheres (DHS). Therein the short-range part of the interaction is 
given by the hard sphere repulsion potential 




(2.17) 




(2.18) 



oo 



for TiA < D 




HC 



= < 



(2.19) 







for Tij > D 



February 2, 2008 



11 



with D being the common hard sphere diameter in the monodisperse case. In the polydis- 
perse case with different hard sphere diameters, Di, one has to replace D by (D,i + Dj)/2 
in the above equation. The only other interaction in this model is the long-range dipolar 
potential. 

Some aspects of the sterical interaction (|2.15jl can also be incorporated into the DHS 
model. Since the magnetic field of a perfect, homogeneously magnetized sphere is exactly 
dipolar one can introduce into the DHS model particles that consist of two concentric spheres: 
an outer hard sphere with diameter D and a magnetized spherical core of diameter D core < D. 
The latter is the source of the dipole field. The nonmagnetic layer models the surfactants, 

n 

i.e., D — D core « 25 po i ymei .. A nonmagnetic dead layer [12| on the particle can also be modeled. 

The DHS model system was initially introduced to describe the electric polarizability of 
fluids of polar molecules. However, in reality deviations from the spherical shape, particle 
polarizability, and higher order moments have an important influence. Ferrofluids, on the 
other hand, resemble much better the model system. That the hard sphere repulsion is 
a good approximation to the Rosensweig potential becomes apparent when looking at the 
quantity exp(— v ) which is more relevant for the equilibrium properties than v itself. In 
Fig.|21exp[— v (r)] is plotted for the Rosensweig potential (solid line) ignoring the irrelevant 
inner attractive region. As can be seen, exp(— v ) resembles quite well the unit step function 
exp(— v HC ) of the hard sphere potential that goes from to 1 at contact distance. 

The long dashed line for the best fitting hard sphere potential and the short dashed line for 
the van der Waals potential are included in Fig. El for comparison. Although the Rosensweig 
potential resembles more the van der Waals potential since both have an attractive part, 
the hard sphere potential fits the Rosensweig potential even a little bit better. 
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III. INITIAL SUSCEPTIBILITY 



Most of the older approaches to the magnetic or electric properties of polar fluids deal 
with the linear response problem, i. e. the initial susceptibility. The nonlinear behavior was 
of lesser interest in common atomic polar fluids since available electromagnetic fields were 
not strong enough to reach saturation. However, with the advent of superparamagnetic 
ferrofluids with large magnetic moments of the particles the linear and also the nonlinear 
magnetization behavior of interacting dipoles became again a more active research area. 

Most of the approaches to calculate the initial susceptibility can roughly be classified into 
cluster expansion, hypernetted chain, and numerical simulation methods. In this section, 
we will mainly address the results for dipolar hard spheres. 



A. Cluster expansion methods 

Consider a system of N particles interacting with each other via a pair potential and 
with an external potential t>j. Then the canonical partition function is given by 

Z = / exp I - Vj - ^ Vjj j dr . (3.1) 

Here dT indicates the integration over the configuration space spanned by the positions and 
magnetic moments of the particles. The ideal-gas factor contained in the canonical partition 
function that comes from the momentum degrees of freedom is not indicated here. 

The key point of the cluster expansion method is to introduce into (|3.1J) the functions 

fij = exp (-Vij) - 1 (3.2) 

which are small when the interaction is weak. Thus, one inserts exp (—■%•) = 1 + fij into 
()3.1j) and expands the integrand in terms of fy. 
Then (|3.1|) reads up to first order for example 

Z = Z° N + Z%_ 2 [ h d«d(j) , (3.3) 
i<j J 

where Z® N is the partition function of the ideal noninteracting iV-particle system. These 
approximations to Z then allow to calculate the thermodynamic quantities of interest. The 
advantage of this expansion is that in contrast to (|3.1|) the approximations to Z such as 
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()3.3|) and higher orders require only low dimensional integrations, that can be performed at 
least numerically. 

In the case of DHS the calculation of the leading terms can be done even analytically 
if one performs another series expansion, namely the expansion of ()3.2j) in terms of the 
dipolar interaction potential v® D that enters via = v® D + vfj C . The two expansions then 
translate into a power series for Z and thus for M or \ in two dimensionless parameters, 
namely the volume fraction 

(3.4) 

of the particles, and a dipolar coupling constant 

2 

TfL 

X = Air^ksT ' (3 ' 5) 

It is also common to choose p* = ND 3 /V and y = 80A/3 as parameters. These parameters 
fully characterize the system for H = 0. To make a comparison with Stockmayer particles 
easier, we will use these parameters here also. In the case of Stockmayer particles one also 
has to specify the potential strength v entering into Eq. (|2.18|) . It should be noted that in 
the literature a system of Stockmayer particles is normally characterized by giving values 
for 1/vq and \/v . 

In ordinary ferrofluids, <ft is of the order of 10~ 2 . For magnetite ferrofluids with 10 nm 
particle diameters and an additional polymer layer of 2 nm A is less than one. In this 
case, higher orders in the (0, A)-expansion have only minor influence on x- But it is also 
possible to produce much stronger interacting ferrofluids. For example, van Ewijk, Vroege, 
and Philipse ^| report the production of ferrofluids with values for A up to 2.7 and volume 
: ractions up to 0.5. Even higher A could be realized by Mamiya, Nakatani, and Furubayashi 

The long-range nature of the dipolar interaction causes some cluster integrals to be ge- 
ometry dependent and mathematically ambiguous if the thermodynamic limit is naively 
performed. This problem is circumvented by dealing with a macroscopic but finite spher- 
ical geometry. In that case the initial susceptibility with respect to the external field 
Xe = dM/ dH e (H e = 0) can be calculated unambiguously and is related to the geometry- 
independent quantity x y i a 

77~ = \xe ■ (3.6) 

3 + x 3 
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In 1966 Jepsen [2Cj performs a cluster expansion to obtain 

x _ i _ A 3 

3 + x 3 XL 144 Xl 
for DHS. x is then up to the calculated order given by 



X Xl •TjXi, 



xl 



(3.7) 



(3.8) 



This equation can also be understood as a (0, A)-expansion, since \l = M sat m/3k B T = 
Nm 2 /3fiQk-QTV can be written as 

Xl = 80A. (3.9) 

The fact that both the Weiss and the Onsager model correctly predict the prefactor of the 
x|-term shows that this term has a mean-field origin being independent of the type of 
short-range interaction. 

Rushbrodce Q crates more *™ ar 19 79 . His result, as corrected later by JosBn 
reads 

160 ,„,„ 8(1187- 600 In 2) 



X 



64 
225 



-0A + — 2 A 4 3 A 3 - -i — ^0 3 A 4 + 1.13358 4 A 3 



3 + x 3 r 225 r 9 3375 

= 2.666670A + O.2844440 2 A 4 - 17.77780 3 A 3 



(3.10) 



1. 827820 s A 4 + 1.13358 4 A a . 



4\3 



The 0A- and 3 A 3 -terms agree with the \l- and x|-terms of Jepsen. The last term re- 
quires a numerical integration. Rushbrookes error affected the 3 A 4 -term. Buckingham and 
Joslin calculate the second dielectric virial coefficient, i. e. the 2 A n -terms in the above 
expansion: 



X 



3 + X 



64 



^3n[(2n + 3)!!]' 



£(3*- 



n 



k=0 



Joslin 



/ \ 

2k 



2 \ 2n+2 



+ 



(3.11) 



221 ] presents besides the aforementioned correction to the result of Rushbrooke also 
the term 



32(315 In 2 - 218) 
3375 ~ 



3 A 5 = 0.003236620^ 



3\5 



Thereafter Rushbrooke and Shrubsall [24] calculate in addition the term 
^105133619 298976 



4042500 



7875 



In 2 



3 A 6 = -0.30* 



3 A e 



(3.12) 



(3.13) 
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A numerical calculation of the 3 -term as a function of A using Monte-Carlo methods is 
reported by Joslin and Goldman in 1993 j^j]. 

Tani et al. take a slightly different road to calculate x using the system of (nonpolar) 
hard spheres as reference system. They get 

64 

X = 80A + y0 2 A 2 + /(0)A 2 , (3.14) 

where the last term is a numerical expression based on approximate expressions for the 
two- and three-particle correlation functions for hard spheres as found in Monte-Carlo 



simulations. For small this result reduces to (|3.8|) . Goldman 27[ compares results of this 
theory very successfully to Monte-Carlo data for systems of hard spheres and Stockmayer 
particles. 

B. Hypernetted chain and related approaches 

Using the hypernetted chain (HNC) approach the initial susceptibility is calculated via 
;he two-particle correlation functions g(l, 2) and h(l, 2) = g(l, 2) — 1 of the polar particles 
. The HNC is based on the Ornstein-Zernike relation 

N f 

h(l, 2) = c(l, 2) + — J c(l, 3)/i(3, 2)d(3) , (3.15) 

defining the direct correlation function c(l, 2) and a closure relation, the HNC approximation 
reading 

c(l, 2) = h(l, 2) - \n[g(l, 2)] - v 12 r 12 > D , (3.16) 

0(1, 2) = r 12 <D (3.17) 



in the case of hard spheres. It can be shown [28j that this approximation is equivalent to 
a cluster expansion using an infinite number of /-integrals, neglecting only a certain class, 



the so-called 
Wertheim 



bridge diagrams that appear only in higher orders. 



29( and Nienhuis and Deutch |30| used the so-called mean spherical approx- 



imation (MSA) [3l| to calculate the initial susceptibility of DHS. In the MSA ()3.16|) is 
replaced by 

c{l,2) = -v 12 for r 12 >D , (3.18) 
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effectively using the approximation ln[g(l, 2)] = ln[l + h(l, 2)] m h(l, 2), which is correct for 
large distances, where h(l, 2) is small. Wertheim showed that this equation together with 
(13. 15 J) allows a simple ansatz, representing the dependence of the correlation functions on 
mi and m 2 as linear combinations of three functions. He sets 



with 



h(l, 2) = h s {r 12 ) + h D (ri 2 )f D {mi, m 2 ) + h&(r 12 )f&{mi, m 2 ) (3.19) 

/ A (mi,m 2 ) = mi-m 2 , (3.20) 
/ D (ihi,m 2 ) = 3(fi 2 • mi)(fi2 • m 2 ) - ihi • m 2 • (3.21) 

(7(1,2) and c(l,2) have a similar representation with the same /d,/a- The MSA then 
provides an implicit solution for x that is given by the three equations 

X = * L ; , Xl = q(2x) - q(-x) , q(x) = ^ + ^ . (3.22) 
q{— x) (1 — x) 4 

Here, \ depends only on \l — 80A but not on <ft and A separately. Expanding \ in powers 
of xl yields again ()3.8(1 for the leading terms. 

Verlet and Weis [22J propose in 1974 an improved theory. Therein they replace the term 
hs(r) that does not depend on the dipolar character of the fluid by a better result for 
nonpolar hard spheres. Stell and Weis 3^] calculate the initial susceptibility of DHS using 
this modification and a further improved version. They get larger values for \ that agree 
better with the Monte-Carlo data. 

Patey improves the theory of Wertheim by retaining the ansatz 1)3.19)1 but expanding 
the full relation ()3.1(ijl linear in fjj and /a- This approach is known as linear HNC (LHNC), 
although it was pointed out later that an equivalent theory was proposed already in 1973 by 
Wertheim j^] himself, called single-superchain theory. In contrast to the MSA, the LHNC 
requires a numerical calculation of the correlation functions hs{r), ho{r), etc. to find the 
initial susceptibility. Within the LHNC \ depends on both <fi and A independently. Patey 
did not apply the "pure" HNC but uses again hs and cs from the system of nonpolar hard 
spheres as system of reference (known as reference LHNC or RLHNC). 

The RLHNC was followed by the reference quadratic HNC (RQHNC) proposed by Patey, 
Levesque and Weis. In j3| this approach is applied to DHS. In the RQHNC (|3.16)1 is 
expanded up to quadratic terms in f D and /a- Differences in \ between RQHNC and 
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RLHNC amount only up to a few percent for A < 2 and relatively large <p ~ 0.4. But they 
are more significant for larger A. However, later on it turned out in Monte-Carlo simulations 
(see Sec. IIII (^1 that both theories overestimate x here — the apparently better RQHNC even 
more so than the RLHNC. 

These authors considered also Stockmayer particles, using DHS as a reference system. 
They argue that a Stockmayer system with = 0.8 ir/6 and l/v = 1.35 should be similar to 
a system of DHS with the same </> for not too large A. Indeed, the results for the susceptibility 
do not differ very much for A < 2. For higher coupling constants the RQHNC now gives 
smaller values for v than the RLHNC. 

I — I 

Agrafonov, Martinov and Sarkisov [3j| expand the full HNC equation in terms of A and 
obtain the result 

X , . A ,2 \2 



^ L = 64^A 2 ^1 - -0AJ (3.23) 

for DHS. This expression yields after solving for x the correct terms in <ft\ and <p 2 \ 2 . However, 
the next term is already incorrect. Chan and Walker propose to expand the HNC 
approximation ()3.16|) given in the form 

g(l, 2) = exp[/i(l, 2) - c(l, 2) - v 12 ] , (3.24) 

in terms of rotational invariants for the orientational distribution of rhi and m 2 , truncating 
the expansion after a sufficient number of terms. They take into account however only three 
terms, namely the functions already used in the MSA, RLHNC, and RQHNC. Using this 
truncated HNC ansatz (THNC) they obtain even higher values for the susceptibility of DHS 
than predicted by the RQHNC. 

In 1985 Fries and Patey [39j finally extend the ansatz (|3.19|) by adding more terms. They 
also use (J3.16|) in a form differentiated with respect to r 12 eliminating the logarithmic term. 
Solving numerically the full RHNC approximation for DHS in that way they come to results 
for x that are in much better agreement with Monte-Carlo calculations. In the same year 
Lee, Fries, and Patey also investigate the Stockmayer system. As in the case of DHS 
the results agree better with numerical data than those of RLHNC and RQHNC. Finally, in 



1986, Fries and Patey 



4l| apply the similar Percus-Yevick approximation to DHS finding 



it less well suited to predict the susceptib ility than the RHNC and even the MSA. 



Another modification is done by Lado 



42j , who varies the diameter of the nonpolar hard 



spheres in the system of reference to minimize the free energy. Lomba, Martin and Lorn- 
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bardero [43] solve the pure HNC finding for the correlation functions even better agreement 
with Monte-Carlo calculations than for the RHNC 



C. Numerical simulations 

Numerical simulations for calculating the susceptibility are mostly based on Monte-Carlo 
methods, x is either determined from the magnetization fluctuations (M 2 ) in the absence 
of an external field or more directly from the magnetization in small fields. Calculating the 
susceptibility from the simulation results for h(l, 2) are less appropriate since x depends on 
the long-range behavior of h(l, 2) that cannot be determined accurately in finite simulation 
cells [36]. 

The influence of the distant dipoles can be incorporated by cutting off the dipolar inter- 
action at a finite distance R s and replacing the dipoles beyond R s by a magnetic continuum 
of given x (reaction field method). This susceptibility is either adjusted to the susceptibility 
of the simulated system or set to an arbitrary value (commonly to % — 0). The chosen value 
then enters into the determination of the system's susceptibility. 

Another method that is the most commonly used today is the Ewald summation based 
on a periodic continuation of the system together with an effective method for calculating 
he resulting total dipolar field. A discussion of the different methods can be found e. g. in 
44lk5j]. The Ewald summation is described in 

Monte-Carlo calculations of the free energy and other thermodynamical quantities that 



the dipolar forces were a 



ready performed 



depend less strongly on the long-range nature o 

in the seventies of the last century j^, 47, 3]- Patey et al. 0, till 5^] calculate the 



correlation function h(l,2) and compare with the predictions of the LHNC and QHNC. An 

indirect comparison for the susceptibility, demonstrating the superiority of the MSA over 

the Onsager model can be found already in 1491 1 . 

An early Monte-Carlo calculation of x is done by Adams and McDonald [46] in 1976. 

The authors however do not consider a fluid but f. c. c. and s. c. lattices. They compare 

Ewald sum and reaction field results for interaction strengths up to A ~ 2. Levesque, 
I I 

Patey, and Weis |5i| calculate x f° r = 0.87r/6 and A < 1. They compare different methods 
and system sizes and in 1982 the y a lso study the case 

A = 2 H- Adams Q considers 
already in 1980 A = 2.75. Lado p2| uses a combination of Monte-Carlo calculation and 
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RHNC for the long-range part of the potential. 

Susceptibility calculations for Stockmayer particles were mainly performed for the values 
<f) = 0.87r/6 and vq = 1/1.35 as already considered by Patey, Levesque and Weis H- Pollock 
and Alder 0] consider values for A up to 3 while Adams and Adams [45] study different 
values for vq. Levesque and Weis calculate \ for A = 2, and Evans and Morriss 0] do 
so for A < 1.7. Neumann, Steinhauser, and Pawley consider the parameter combination 
Vq = 1/1.15, 4> ~ 0.43, and A = 2.6. Hesse-Bezot, Bossis, and Brot |5jj use in 1984 molecular 
dynamic methods to investigate A = 2. 

Some of the works cited above treat the long-range nature of the dipolar interaction 
i—,,, see the di_ by N_ 0, Ne Umam and Steinnanser 0Q. and 
Gray et al. 0] • These papers also discuss the relationship between the reaction field and the 
Ewald sum approach. Further discussions of the Ewald summation and of other techniques 
used for dipolar systems were given recently by Wang, Holm, and Miiller 0, 0] . 



D. A comparison 

Unfortunately most of the older works deal with systems with large volume fractions: a 
typical value is <fi = 0.87r/6 = 0.419. This is reasonable for customary polar liquids but it is 
very high for ferrofluids. Fig. H] shows a comparison of the most important theories discussed 
above for the case of DHS. In this figure % is plotted as a function of A for fixed (f) = 0.419. 
The Monte-Carlo data (denoted by circles) and the results for RLHNC and RHNC were 
taken from [39] and a similar plot in 0]. Onsager theory and MSA clearly underestimate 



the Monte-Carlo data while Weiss theory and RLHNC overestimate them. The Rr 



NC 



• 126]. 



works much better, as does the numerically much less complicated result by Tani et al. 
The curve denoted as "Huke and Liicke" refers to a cluster expansion theory |£J] for the full 
magnetization curve. In the linear case considered here in this section this theory yields 
Eq. ()3.8|) and in addition also the <j) 2 \ n terms given by Eq. ()3.11[) . 

The large scatter of the Monte-Carlo results for A = 2.75 is indicated in Fig. [3] as well. 
Surprisingly there seems to be no final conclusion in the literature concerning the appropriate 
value of x f° r ^ = 2.75 despite the fact that these different Monte-Carlo results have been 
in the literature for quite a while. The different predictions for % at A = 2.75 are discussed 
6f| where preference is given to a relatively low value of x ~ 64. 



m 
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E. Other potentials 

The susceptibility of some variants of the DHS and Stockmayer particle systems were also 
investigated. The system of dipolar sticky hard spheres, where the short-range potential 
contains an additional 5-function term, were investigated by Chapela and Martina |fi5j using 
integral theories like MSA and RLHNC to calculate x- The authors come to contradictory 
results concerning the effect of the additional term in the potential on \. Joslin and Gray 
jf^l show that the second dielectric virial coefficient of dipolar sticky hard spheres is positive. 

Kusalik investigates in 1989 and the following years dipolar soft spheres, i. e. 

Stockmayer particles without the attractive r _6 -term in the short-range part of the potential 
using the RLHNC and the RHNC as well as Monte-Carlo calculations. The Monte-Carlo 
results for v = 1/1.35, = 0.87r/6, and A = 2 are very similar to those for Stockmayer 
particles. But RLHNC and RHNC overestimate x somewhat more. 

Henderson, Boda, Szalai, and Chan jzQ, 71] apply the MSA, the theory by Tani et al. 0], 



and Monte-Carlo calculations to dipolar Yukawa particles with a hard core. Again, the 
results are found to be similar to those for dipolar hard spheres and Stockmayer particles. 
The Monte-Carlo data for x at </> = 0.87r/6 and A = 1, 2 depend little on the strength of the 
Yukawa potential and they agree well with the Monte-Carlo data for dipolar hard spheres. 
The perturbation theory by Tani et al. is found to reproduce the Monte-Carlo data better 
than the MSA. 

To summarize, the nature of the short-range interaction seems to have a rather weak 
influence on the susceptibility, even for relatively high densities. 



F. Polydisperse theories 

In real ferroffuids the influence of polydispersity can normally not be neglected. In the 
case of ideal paramagnetism discussed in Sec. Ill Al the initial susceptibility is given by 

^=Es^b' < 3 - 25 > 

as follows from (J2.5|) . Therefore xl is proportional to the second moment of the distribution 
of the magnetic moments and larger than in the monodisperse case if the average particle 
volume is kept constant. When particle interactions are taken into account even higher 
moments of the distribution enter into the initial susceptibility. 
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Several of the theories discussed above were also extended to the system of polydisperse 
or at least bidisperse DHS, where the particles have different hard sphere diameters and/or 
carry different magnetic moments. Already in 1973 the MSA of Wertheim is applied to 
the polydisperse case by Adelman and Deutch although restricted to particles with a 
common diameter. The authors define an equivalent monodisperse system having the same 
thermodynamical properties. 

Isbister and Bearman [x^] generalize the MSA to mixtures with arbitrary diameters. The 
susceptibility predicted by this theory is calculated by Freasier, Hamer, and Isbister j^J. 
The authors also describe a comparable monodisperse system gi ving at least approximately 
the same results. Ramshaw and Hamer simplify the evaluation [75]. Cummings and Blum 
761 compare results of the MSA to Monte-Carlo results for the bidisperse case. 

n n 

Lee and Ladanyi present |7j| a RLHNC extension and generalizations |7j| of the RHNC 
and the cluster expansion by Tani et al. j^J. They compare to the Monte-Carlo data 
from [2]], coming to similar conclusions about the quality of the different theories as in the 



monodisperse case discussed above (Sec. IHI Dj) . 
G. Applications to ferrofluids 

With the preparation of highly concentrated ferrofluids real physical systems became 
available that have much more in common with the simple model system of DHS than the 
ordinary polar fluids. Ferrofluids thus allowed not only more reliable experimental tests but 
they also initiated new theoretical investigations. 

Two-dimensional Monte-Carlo simulations on the initial susceptibility of ferrofluids are 
performed by Menear, O' Grady, Bradbury et al. in the early eighties Here 
the particles are not modeled as pure hard spheres but the sterical repulsion of the surfactants 
according to Rosensweig (|2.15jl is also taken into account. In Q| a polydisperse ferrofluid 
is simulated. In particular, the temperature dependence of the susceptibility entering via 
A ~ 1/T is investigated. A Curie- Weiss law I fx ~ T — T is found as to be expected from 
simple mean field models, e. g. the Weiss model. 

Bradbury, Martin, and Chantrell also perform full three-dimensional simulations 

for DHS with additional sterical repulsion and obtain good agreement with RHNC calcula- 
tions. The authors again find a Curie-Weiss law for x{T) but less pronounced than for the 
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2D calculations. The calculations are done for a small = 0.01 and very high A up to 7.21. 
These values refer to the real hard sphere diameter. However, the dipolar coupling would be 
better characterized by using the total diameter of the particles, i.e., of core plus polymer 
layer. That would reduce A to values of maximal 2.63. In this A-range much simpler cluster 
expansion theories still work and reproduce the result. It should be mentioned that although 
1/x — c(T — To) is a very good approximation in that range, the Weiss model already fails 
since it predicts significantly different values for c and T . 



79 



Q- The results in Q 



A behavior of 1 jx ~ T — T is also found in experiments 
however rely on data in a narrow range of temperatures of 60 K where only very strongly 
interacting ferrofluids should exhibit a pronounced nonlinear behavior. In |^| a negative 
value for To is found, i. e., a Neel behavior that was not predicted by the theories. A more 
comprehensive experimental investigation of the initial susceptibility is published in 1990 
by Holmes, O'Grady, and Popplewell 87]. Inspecting a broader range of temperatures the 
authors show deviations from the Curie-Weiss law for ferrofluids based on magnetite. They 

n 

also argue that the Neel behavior found in [8Qj may result from measurements for too large 



fields, i. e. already in the nonlinear range of the magnetization curve. They find a linear 
and a quadratic term in the concentration dependence of the susceptibility. The nonlinear 
behavior of l/x^lis confirmed by Williams et al. Ssl in 1993. 
Morozov et al. 



89j,|90j apply the MSA to ferrofluids. Further results can be found in [91 1 



The authors demonstrate good agreement with experiments concerning the concentration 



and temperature dependence of x- As discussed above in Sec. IIII D| the MSA does not seem 
to be a good theory for strongly interacting systems. The investigated ferrofluids {(f) < 0.172, 
< 10) however fall in range where the difference in x between MSA and cluster expansions 
is only a few percent. Applying the formula ()3.8|) shows that the x|-term which is 
correctly predicte d by the MSA gives only minor contributions to x- 

Pshenichnikov |9_3J] compares several theories, namely the Onsager model, the Weiss 
model, the MSA and a theory by Ivanov et al. |94] with experimental data in a similar 
range of and A. The latter theory deals with the equilibrium magnetization for ar- 
bitrary fields. It is explained in the next section. For the initial susceptibility this theory 
gives x — Xl + x|/3; l - e -> Eq. (j3.8j) up to second order. Pshenichnikov shows good agree- 
ment with both Ivanov's theory and the MSA, thereby also indirectly demonstrating the 
small influence of the x|-term that is present in the MSA but not in Ivanov's theory. Both 
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theories however fail to adequately describe the temperature dependence of ferroffuids with 

x{T) lying in the range between « 20 and 80. 
I I 

Kalikmanov [9_2j proposes in 1999 the so-called algebraic perturbation theory that is 
corrected by Szalai, Chan and Henderson [3| and found to be identical to the cluster 
expansion theory of Tani et al. 



Van Ewijk, Vroege, and Philipse 18] perform susceptibility measurements on highly 
concentrated ferroffuids and compare with different theories. For a ferroffuid with A ~ 2.7 
and xl < 6, they find a surprisingly low susceptibility of x < 10, best described by the 
Onsager theory or the MSA. These results seem to be in contradiction to the data presented 
in Fig. HI although <p is slightly higher there. Some theoretical support to these findings can 
however be found in a paper by Pshenichnikov and Mekhonoshin who present in 2000 
Monte Carlo data in the same range of xl- For A = 3 and A = 4, they find similar low values 
for the susceptibility. They explain the low susceptibility with aggregates of particles, that 
have a small combined magnetic moment and interact with the field only weakly. 
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IV. MAGNETIZATION 

The equilibrium magnetization of a ferrofluid cannot be described in general by a sim- 
ple Langevin ansatz or its polydisperse generalization. Experimental hints to this fact can 
already be found in a paper from 1979 by Tari et al. 0| where a comparison of a magnetiza- 
tion curve with a theoretical curve based on an independently measured size distribution is 
made. However, the theories for the magnetization of, say, dipolar hard spheres in arbitrary 
fields are less advanced than for the special case of the initial susceptibility. Most of the 
theories discussed below have been devised with regard to the application to ferrofluids. For 
this reason, many of them were directly or shortly after extended to include also the effects 
of polydispersity. 

When comparing a polydisperse ferrofluid with a monodisperse one with the same av- 
erage magnetic moment of the particles and the same volume fraction, the equilibrium 
magnetization curve will be steeper for small fields in the polydisperse case since the initial 
susceptibility is higher. The asymptotic behavior for high fields will however be the same. 
This results in a maximal difference between the magnetization curves for medium fields 
a« 1. 

Pollock and Alder [5^ and Adams and Adams perform Monte-Carlo calculations for 
the equilibrium magnetization of Stockmayer particles. In 1981 H0ye and Stell 3| propose a 
mean field model based on general thermodynamic considerations and compare their results 
with the data from [53]. The magnetization in this model is given by 

M = M sat C ( -^-H eS ) , (4.1) 



where H efi is an effective field 



H cS = H+^(l-Q)M . (4.2) 



Here is to be taken from the results for the initial susceptibility. For O = the formula 
reduces to that of the Weiss theory. Dikanskii [99] proposes a concentration depen dent 
term for an effective field calculated from a fit to susceptibility data. Sano and Doi jlOo| 
investigate phase transitions in ferrofluids using a model of a randomly filled cubic lattice. 



As an additional result they o 
Woodward and Nordholm 



jtain the Weiss expression for the equilibrium magnetization. 



lOlj l propose an effective potential for the dipolar interaction 
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that results from averaging over the orientations of the dipoles according to 

exp (-Vg/ksT) = <exp (-V$ D faT))^ . (4.3) 

The so defined V is then a functio n of p article separation and temperature. It vanishes 
for large distances as rf 2 6 . In 1986 [102] these authors extend their ansatz to a theory 
for the equilibrium magnetization. They propose a functional for the free energy F of 
a magnetic continuum in a magnetic field. Minimized with respect to the orientational 
distribution of the dipoles it gives for zero field an expression for F that contains V cS as 
interaction potential. Including the external field the dipolar interaction is described by a 
long- range part as in the Weiss model and a short-range part that is a generalization of 
V eS which, however, depends also on the magnetic field. The magnetization has the form 
()4.1I4.2|) with a function B(M, T) that can be calculated numerically. The authors compare 
the resulting initial susceptibility with the RLHNC predictions. The results are however 
not convincing: They are larger than in the RLHNC already for <fi m 0.3 and diverge for 
<p 0.4. Com parisons for DHS with the magnetization results of H0ye and Stell can be 
found in [mJ. 

Berkovsky, Kalikmanov, and Filinov I104I llOsj ] develop between 1985 and 1992 a 
thermodynamical theory for ferrofluids. Based on cluster expansion methods they calculate 
the magnetization in two special cases. They derive expressions for spherical geometries 
where H = H e — M/3. For small fields a < 1 the authors obtain for the susceptibility with 
respect to the external field H e 

dM , , 160,o ,0 , 
Xe = qjj- = 8A0 - — A 3 3 . (4.4) 

If one restates this expression into a result for \ a missing 2 A 2 -term appears and one 
obtains (|3.8|) . For stronger fields they get 

M = M sat [C(a e ) + L li2 (a e )0A 2 /(0)] . (4.5) 

Here Li^{a e ) is an analytic expression and /(</>) depends on the correlation function g(l, 2) 
of nonpolar hard spheres that is used as reference system. Again, by restating the expression 
in a form M = M(H) an additional term of order <fr\ shows up. Here one has to note that 
M/M sat is independent of and A in leading order, whereas \ starts with a O(0A)-term. 
Likewise, terms of order m A n in M/M sat always correspond to O(0 m+1 A n+1 )-terms in \. 
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. 1Q7, 108,0 



Between 1990 and 1992 also Ivanov et al. [94J, llM llOTj, UOgJ, [109| develop a thermody- 
namical theory of ferrofluids also using cluster expansion techniques. For a needle-shaped 
geometry in which H e = H they obtain the correct result in 0(<j)\) 

M = M sat [C(a) + 84>\C(a)C'(a)} . (4.6) 

A correct polydisperse generalization, replacing £(a) by £ poly (a) ()2.7|) . is also given. 

Lebedev |ll0j suggests to apply the MSA also to finite magnetic fields. To that end he 
replaces the equations ()3.22|) by 

M L (H) dM(H) q(2x)-q{ 



M(H) 



-x) 



(4.7) 



q(— x) ' dH q(—x) 

with q(x) given as before by the third equation in ()3.22j) . Thus, \ and xl in the first 
equation of ()3.22|) is replaced in ()4.7|) by M and M L , respectively. The second equation of 
()4.7|) can be derived by combining the first two equations in (j3.22j) and replacing the initial 
susceptibility x by the susceptibility dM(H)/dH for finite fields. Hence, the equations ()4.7|) 
reduce to the MSA for sma ll fields. 

Morozov et al. 0, llll| present in 1990 a real extension of the MSA aimed at predict- 
ing the magnetizat ion for arbitrary fields. They use the so-called Lovett-Mou-Buff-Gubbins 



equation 



112 



113| to relate the one-particle probability density p(r, m) in the presence of an 



external field to the direct correlation function c(l, 2) taken in the MSA form of Eq. (|3.18|) . 
With Wertheim's method 2^| they then derive from the Ornstein-Zernike equation two inde- 
pendent Percus-Yevick like equations for a hard-sphere fluid with renormalized densities. As 
a result of this approximation scheme they present M given by the Langevin magnetization 
for an effective field Ca 

M = M sat C (Ca) 

Here C has to be calculated as follows: Let 
<91nM 



(4.8) 



.4 



and 35a; 2 -5(10 -3A)x - 7(7 + 3 A) = 



(4.9) 



d\n(CH) 

With xi and x 2 being two solutions of the quadratic equations, C must be given selfconsis 
tently by 



C 



1 + T (1 + Xi) h\ q{ll y) 



4vr . .Ml 1 , 
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Here y is defined via the second equality and q(x) is still given by the expression in (|3.22|) . 
The authors also propo se a polydispersive generalization. Comparisons with experiments 
can be found in |9lLlllll |: The MSA reproduces well the magnetization curves of polydisperse 
ferrofluids with saturation magnetizations up to 87 kA/m. 

Based on preliminary wor k yJJf concerning the 2D case Bradbury and Martin |115| and 
independently Ayoub et al. jll6j | develop the so-called dimer model of dipolar interaction. 
The partition function for N particles is written as the (iV/2)th power of a two-particle 
partition function calculated via a cluster expansion assuming that the two particles are 
closer than a typical next-neighbor distance. This model has several shortcomings. The 
long-range character of the forces is not taken into account leading to an expression for M 
in which the 0A-term is missing. Also the calculated O(0A 2 )-terms are wrong in both cases 
and moreover they do not agree with each other |m|. Even the corrected calculation leads 
to a wrong result by a constant factor due to the somewhat arbitrary assumption of using 
the next-neighbor distan ce as a cutoff. 

Zubarev and Iskakova jl present a magnetization equation relying on a model of non- 
interacting magnetic clusters (chains). They derive the equilibrium distribution of chain 
lengths based on a simple expression for the free energy of a system of such chains. The 
magnetization is then formally the equilibrium magnetization of a polydisperse ideal param- 
agnetic gas where the chains are the fundamental constituents. The chain length distribution 
is howeve r a f unction of the magnetic field and furthermore of and A. Abu-Aljarayesh and 
Migdadi |l lo| investigate this theory further, calculating the entropy and some equilibrium 
magnetization curves. 

Pshenichnikov, Mekhonoshin and Lebedev jl2o| | suggest an equation for the equilibrium 



magnetization that reads 



M = M sat £ 



mM sat „, 

ot H : — — £(a) 

3k B T K ' 



(4.11) 



It agrees better with their experiments in the range of medium fields than the expression 
(14. 6(1 that gives the same initial susceptibility. It furthermore gives more consistent results 
for the diameter distribution of polydisperse ferrofluids at different stages of dilution. Fur- 
thermore, their earlier result (|4.fij) can be derived from Eq. (I4.11j) by expanding the latter 
up to first order in <j) and A since mM sat /3k-QT = 80A. On the other hand, Eq. (J4.ll)) can be 
obtained from the Weiss formula (|2.8|) by iterating and stopping after the first step. How- 
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ever, Eq. (J4.ll)) has the advantage over Eq. ()4.6|) that it always produces physically sound 
magnetization curves with M being a monotonic function of H with < M < M sat . This 
is not the case for Eq. (J4.6j) when <p\ becomes too large. 

In 2000 the authors of this review 3] perform a cluster expansion calculation to derive 
some of the terms in the sum 



M — M, 



sett 



£(a) + ^L a , fe (a)0 a A fe 



(4.12) 



namely a number of L^-terms and the leading term, L 2j 2, in O(0 2 ). All the calculated 
terms are given as analytic expressions. The first one, Li i i(a) = 8£(a)£'(a), yields the 
result ()4.6|) of Ivanov et al. and L\ ? is th e same as in the expression ()4.5|) . A generalization 
to the polydisperse case is given in |l2l| . 

In the same year, Klapp and Patey |l22| present results obtained from two different 
integral equation approximation schemes that use either the MSA or the RHNC as two- 
particle closures. The singlet density p and the external field is related in formally exact 
ways either via the first BGY hierarchy equations |28| to th e pair correlation function g(l, 2) 
or via the Lovett-Mou-Buff-Wertheim equations |l!2 . Il2 j ] to the direct correlation function 
c(l, 2). The results turn out to be strongly dependent on the form of the closure used. For 
example, the authors find in zero fiel d ferr omagnetic behavior using the RHNC closure that 
is absent in case of the MSA closure 122]. 

Pshenichnikov and Mekhonoshin [96] present Monte-Carlo results for a finite system 
giving x as a function of xl for fixed A = 1, 3, and 4. They also calculate the equilibrium 
magnetization for a volume fraction = 0.226 and A = 1 and 4. The results for A = 1 
can quite well described by the equation 1)4.11)1 . The equilibrium magnetization for A = 4 is 
however smaller than predicted by any theory for interacting particles and is explained as 



result of the formation of c 
Ivanov and Kuznetsova 



usters with a small total magnetic moment. 



124J perform a cluster expansion calculating the O(0 2 Af]__term 



in M/M S n.t.. While their resu lts a re not given explicitly they seem to reproduce [12q the 
result of 0] . The authors of |l24| suggest an extension of Eq. 1)4.11)) that reads 
'mH cS " 



M = M sat £ 



(4.13) 



k B T J ' 3 ° v ' 144 sat v dH 

A polydisperse generalization, replacing C by £ po i y is also given. The formula ()4.13)) produces 
sound magnetization curves for arbitrary <f), A as in the case of Eq. ()4.11j) . When comparing 
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the results of 



1241 ] with those of the cluster expansion 



64fl one has to be aware that in 



124 1 



some calculated terms that do not contribute to the initial susceptibility were discarded in 
order to be able to write the result in the form (J4.13)) . Thus, when this equation is expanded 
in 0A one does not recover the general O(0 2 A 2 )-term of 1)4.12)1 while it is reproduced in the 
expression for the initial susceptibility. 

Molecular dynamics simul ation s of monodisperse dipolar systems are performed in 2002 



by Wang, Holm, and Miiller 



12q for several combinations of <p and A. They use a modified 



Stockmayer potential with a shifted and tr unca ted Lennard- Jones potential. Comparing 
with the theory of Ivanov and Kuznetsova |l24| they find good agreement for A < 2 but 
deviations for higher A, both for the magnetization curve and the initial susceptibility. 
The authors also demonstrate that M and y are not functions of the product 0A alone as 



predicted in 



124j and that the theory 124j becomes less applicable for large A and small 



0. Furthermore they analyze the microstructure of the simulated fluid and determine the 
length distribution of chains. They show that the susceptibility can be better explained 
by replacing the Langevin susceptibility in (J3.8)) by the susceptibility of a system of ideal 
dipolar chains with a length distribution found in the simulations . In 2004 Ivanov, Wang, 



and Holm 



127j use a variant of the chain formation model of jll8^ to obtain a length 



distribution independent of the simulation results. Alre ady in 2003 Wang and Holm publish 
simulation results on bidisperse dipolar systems \l2^ . For a system containing particles 
with A = 1.3 and A = 5.32 and a to tal volume fraction of <fi = 0.07 they report that the 
theory of Ivanov and Kuznetsova |l24l | works well for volume fractions (pi < 0.02 of the large 
particles. However, differences become apparent for larger fractions of the large particles, 
say, 4>l = 0.05 and even more so in the limiting case, <Pl = 0.07, without any s mall ones, i.e., 



for a monodisperse fluid containing only large particles. Kristof and Szalai jl29j | perform 
in the same year Monte-Carlo simulations for two polydisperse ferrofluids with realistically 
chosen paramet ers. T hey find good agreement with the polydisperse version of fj4.6jl . 

Ra§a et al. measure in 2002 the initial susceptibility and the magnetization of 

ferrofluids at various stages of dilu tion. The susceptibility is found to be best described by 
Ivanov's theory |94j and the MSA 



ml 
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V. NONSPHERICAL POTENTIALS AND HIGHER-ORDER MAGNETIC MO- 
MENTS 

Real ferrofluid particles can be expected to have small deviations from a perfectly spheri- 
cal shape. Such deviations affect not only the form of the nonmagnetic part of the interpar- 
ticle potential. Also the magnetic interaction potential will be changed since the magnetic 
field of a nonspherical magnetized particle will contain also higher magnetic moments be- 
yond the dipolar one. Both features will have some influence on the susceptibility and the 
magnetization. The susceptibility of systems of simple nonspherical dipolar particles and 
particles with dipolar and quadrupolar particles where also investigated on some occasions. 

A straightforward gener alizat ion of the system of DHS is the system of dipolar hard 



ellipsoids. Joslin calculates 1311 ] numerically the second dielectric virial coefficient for this 
system for different A and axis ratios. For spheres, this coefficient is positive, see Eq. (j3.11|) . 
However, it depends strongly on the axis ratio: it grows with growing prolateness but it 
becomes negative already for only slightly oblate particles. Perera and Patey [lSlJ ] investigate 
this system by means of the HNC, finding the same trend. 

Dimer ferrofluid particles were modeled by dipolar fused spheres (dumbbells). A Monte- 
Carlo simulation of Stockmayer dimers can be found in a paper of de Leeuw and Quirke 133| . 
X is found to decrease with growing sphere sepa ratio n. Dimers of dipolar hard spheres are 
simulated by Lomba, Lombardero, and Abascal [134]. They study the density dependence 
of the susceptibility for fixed dipolar strength and find a decrease of \ at high densities due 
to the hindrance of particle rotation. 

Dipolar particles of spherical shape that carry, however, an additional quadrupolar mo- 
ment are investigated with respect to their susceptibility by Patey, Levesque, and Weis jl35| | 
and Carnie and Pa t e y Q for differ choices of th e q nad rUpol a r tensor using integral 
theories. The susceptibility is shown to decrease with growing quadrupolar moment. 

In view of the fact that the average shape of a ferrofluid particle is roughly spherical it 
seems to us that systems of identical nonspherical or multipolar particles that have been 
considered in the literature are probably less applicable to describe real ferrofluids than 
systems of DHS or Stockmayer particles. To improve upon monodisperse models of identical 
particles one would have to consider at least binary mixtures of particles of different nature, 
e. g. of prolate and oblate ellipsoids or of monomers and dimers. 
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VI. PHASE TRANSITIONS 



Two different kinds of phase transitions are addressed in this section. The first one is 
the condensation of magnetic particles into a denser phase. The second one is the transition 
from an isotropic nonmagnetic phase to a state with spontaneous long-range magnetic order. 



There are other, less intensive 
such as fl uid / solid transitions 



systems 



137 



y investig ated transitions in systems of dipolar particles 



138 



139j | or the demixing of particle species in bidisperse 



140J. Another example of a phase transition that we do not touch here is the 



emergence of a dipole glass state at lo w temperatures for systems of dipoles in a solid matrix, 



like frozen ferrofluids |14lL 1421 . Il43| . The existence of such a state in liquid ferrofluids has 
been reported, e. g., in jl44 . But then it was shown that the observed temperature behavior 
of the susceptibility can be e xpla ined as an effect of the temperature dependence of the 
viscosity of th e carrier liquid 1451 ] . A review of this subject for the electrical case can be 



found in 



146|. 



A. Dilute/dense phase transitions 

Ferrofluid particles possess an isotropic attractive potential: the outer tail of the van der 
Waals attraction that is not shielded by the steric or electrostatic repulsion. This can result 

n 

in the creation of small clusters, as seen by Donselaar et al. in cryo-TEM images 147 . 

When the attractive potential has a strength of about k^T a dilute/dense phase transition 
should occur. It is possible to weaken the repulsive shielding of the ferrofluid particles enough 
to allow such a reversible phase transition without reducing it so much that an i r reve r sible 



148, 



149 



150 



aggl omeration takes place as observed in experiments with ionic ferrofluids 
[l5jj . This reversible transition is a kind of gas/liquid transition. It affects however only the 
magnetic subsystem of the ferrofluid - the carrier system remains always liquid. We will use 
the term "dilute/dense" instead of "gas/liquid" phase transition also when referring to the 
common model systems without any carrier liquid to avoid a confusion of terms. The dense 
phase exists in the form of (mesoscopic) droplets. 

In a typical ferrofluid the isotropic interaction is rather weak. In Fig. |3J for example, the 
best fitting Lennard- Jones potential has a depth of vq = 0.16. S ince a dilute/dense transition 



does not occur in pure Lennard- Jones systems when vq < 0.76 |152| one concludes that such 
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a phase transition might possibly occur in ferrofTuids with smaller vq only if the dipolar 
interactions have a net supporting, i.e., attractive effect. 

A dilute/dense transition observed in experiment that might be caused mainly by dipo- 
lar interactions was reported by Mamiya, Nakatani, and Furubayashi jlfll ]. These authors 
investigate a ferrofluid consisting of iron-nitride in kerosene with 3 ^$ A ^$ 10. They find 
a hysteresis in \ as a function of temperature when is roughly between 0.002 and 0.06. 
For larger (f> the hysteresis vanishes. Furthermore, x{T) is larger when cooled than when 
reheated in the hysteresis region. The authors explain this with the presence of closed rings 
or drops of particles with mostly vanishing total magnetic moment in the low-temperature 
phase. 

When studying phase transitions of dipolar hard spheres, the dipolar coupling constant 
A ~ 1/T serves as the dimensionless inverse temperature. Likewise, serves as dimensionless 
density. In the case of Stockmayer particles both, A and v , provide an inverse temperature 
scale and in general the latter is used. The other quantity or the temperature independent 
quantity X/vq serves as an additional parameter. 

Considering these main model potentials it is clear that a system of Stockmayer particles 
will exhibit a dilute/dense phase transition below some critical T c since it contains the 
Lennard- Jones potential as special case. The additional dipole-dipole interaction raises T c 
(see below) an effect that cannot be explained trivially, since this interaction can be both 
repulsive and attractive and its orientational mean is zero. However, an averaging procedure 
as in ()4.3|) that takes into account the fact that attracting configurations are preferred yields 
an attracting mean dipolar contribution. 

In the case of dipolar hard spheres and also soft spheres it is even more unclear a priori 
whether such systems show a dilute/dense phase transition since here the dipolar interaction 
is the only one that could cause such a behavior. As we will discuss below, the question 
whether this phase transition exists or not in pure dipolar hard sphere systems is still being 
debated. 

In systems of Stockmayer particles with small vq strong dipolar couplings seem to be nec- 



essary in order to raise the critica 
For example, Frodl and Dietrich 



temperature enough to observe a dilute/dense transition. 
I53I ] find using a density-functional theory that already a 



dipolar coupling constant of A ~ 1.5 is necessary to lower the critical point to about Vn ~ 0.4 



According to Monte-Carlo calculations by Hendriks, van Leeuwen et al. 



154. 



155 



be- 
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tween 1989 and 1994 it is A ~ 2.35 for vq ~ 0.4, and Monte-Carlo calculations of Stevens and 
Grest 157| from 1995 give a very similar value here. That means that common ferrofluid 
particles with small vq will not show a dilute/dense phase transition within a reasonable 
temperature range except in systems where the dipolar interaction is the by far dominating 
contribution. It may seem that in such a case the short-range attraction should play no 
more any role and that Stockmayer particles with high A should behave similar to dipolar 
hard spheres. However, as mentioned already this does not seem to be the case. 

In the 1970s and 80s the dilute/dense phase transition of dipolar hard spheres was studied 
by using mainly analytical and semianalytical theories that were conceived to investigate 
the equilibrium magnetization or initial susceptibility. A usual kind of (V, T)-phase diagram 
was found with a roughly parabolic coexistence curve and no phase transition above some 
critical T c (i. e. below some critical A c ). However, the calculated critical values differed 
somewhat. 



Rushbrooke, Stell, and H0ye 158] present in 1973 the coexistence cur ve ba sed on the 
MSA and a perturbation expansion ansatz by Stell, Rasaiah, and Narang 15^ ] . An early 



Monte-Carlo calculation is undertaken by N g, Va lleau, and Torrie in 1979 160 



N = 32 particles. Woodward and Nordholm 



using only 



10 lj use three similar expressions for the free 



energy ba sed o n their effective potential (14.31) to calculate the coexistence curve. Kalikmanov 

nn 

et al. Ill M05J use the same effective potential and replace it by a related Lennard- Jones 
potential to find the critical point. Joslin and Goldman [^J finally calculate in 1993 the crit- 
ical point by using second and third virial coefficients obtained from numerically calculated 
cluster integrals. The approximate results for the critical point (0 C , A c ) are summarized in 
the following table: 



Work 


[158] (MSA) 


[158] (SRN) 


[1601 


[101]a 


[101]b 


[101]c 


[13] 


[25] 




0.056 


0.083 


0.15 


0.17 


0.14 


0.13 


0.16 


0.028 


A c 


4.4 


3.6 


4.0 


2.3 


2.9 


2.9 


1.5 


3.8 



Despite the different predictions for the critical point there was a general belief at that 
time in the existence of a first-order dilute/dense phase transition. This changed when a 
couple of new Monte-Carlo simulations failed to detect such a transition. 
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Caillol 



16 1| performs a Monte-Carlo simulation with up to 512 particles and does not 



find a phase transition for two coupling constants A = 4.5 and A = 5.55 and volume fractions 



137 



162] 



0.08| < < 0.38-. Furthermore, Monte-Carlo calculations by Weis and Levesque 
done for systems of similar size do not show a dilute/dense transition at low densities, 
0.01 < <p < 0.16, even for couplings as strong as A = 12.25. Instead association of particles 
into chains is observed. 

Van Leeuwen and Smit jl63| consider a dipolar fluid with a modified Lennard- Jones 
short-range potential of the form 



SR 



4wn 



D 



12 



( D 



- I -c - 



(6.1) 



The introduction of the parameter c allo ws a smooth transition from dipolar soft spheres 
(c = 0) to Stockmayer particles (c = 1) For one fixed vq the authors do not find a 

dense-dilute phase transition if c < 0.3 and conclude that the dipolar soft and hard sphere 
fluid will not exhibit s uch a transition. They observe chain formation instead. 

Levesque and Weis jl65| investigate the system of DHS with a dipolar coupling A = 12.25 
in an even larger range of volume fractions, 0.005| < < 0.8|, and in addition the range 
4 < A < 12.25 for = O.Olf and = 0.3f . For A = 12.25 they find long chains when 
< 0.2|. For higher volume fractions the situation is less clear. Shorter chains can still 
be found applying an energetical criterion, but these chains are not directly visible in the 
configuration snapshots. At = 0.6 1 and higher they find a ferroelectric state (see below). 
They observe a continuous association of monomers into chains when increasing A for fixed 
volume fraction. 

In 1994 and 1995 Stevens and Grest 157 . Il66| perform Monte-Carlo calculations for 



dipolar soft spheres both with and without applied magnetic field. They do not find a phase 
transition for the zero field case. Instead they observe also in t his sy stem chain formation 



similar to that of DHS. However, still in 1995 the same authors 



152j | report a dilute/dense 



phase transition for Monte-Carlo simulated Stockmayer fluids where the critical temperature 
measured in 1/vo and the critical density depend linearly on A/% The cri tical parameters 



154, 



i3, r 



155, 156] who considered 



agree well with those found by Hendriks, van Leeuwen et al. 
a smaller range of X/vo already between 1989 and 1994. 

The association of ferrofluid particles into chains for A S> 1 and low densities is a natural 



consequence of the highly directional dipolar potential 



167, 168]. The magnetic moments 
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should be oriented in these chains in the energetically most favorable head-to-tail con- 
figuration with energy ~2Xk^T for a pair of them. The numerical observation of chain 
formation instead of a homogeneous dilute/dense transition leads to a number of papers 
considering simple thermodynamic models of chai ns of dipolar hard spheres of similar parti- 
cles. Small-angle neutron scattering experiments |l6ij ] have provided experimental support 
for the existence of such chain structures. Recently, a direct visual confirmat ion of the 
existenc e of chain structures in cryo-TEM images was achieved by Butter et al. jl7o| |. 

Sear jl71 | argues that the chains should interact with each other only weakly and so he 
considers a model of only nearest-neighbor inte ractin g pa rticles, for ming ideal chains. He 



compares with results of the simulations of jl65j | and 16 1| . Van Roij |l72| investigates the 
competition between condensation due to isotropic interactions and chain formation due to 
dipolar forces. 



The most detailed model was proposed by Osipov, Teixeira, and Telo da Gama 



173. 



who also take into account chain-chain interaction. They compare with results from 161 | 



and |162| . These autho rs find semi-quantitative agreement concerning the critical points 
with the results of jl63| for Stockmayer-like potentials. They predict a dilute/dense phase 
transition to appear also for particles with weak isotropic attraction (but not for DHS) where 
in both phases chaining should be observed. However, they also note that finite-size effects 
may make it difficul t to o bserve this kind of transition in simulations. In 1999 Tavares, Weis, 
and Telo da Gama 17^ improve the theory and perform own Mo nte-Carlo calculations at 
low densities = 0.05| and 5 < A < 7.6. Also in 1999 Levin jl76| uses the Debye-Hiickel- 
Bjerrum theory to show that the clustering of particles into chains in the system of DHS 
should always reduce the density of free particles below the necessary density for a phase 
transition. 

Klapp and Forstmann 



177 



178j investigate Stockmayer and DHS fluids using the RHNC. 
There exists a parameter region of low temperatures where the RHNC fails to provide 
solutions. The boundary of this region is believed to be connected with phase transitions. 
At small densities the authors find growing fluctuations near this boundary that seem to 
represent the appearance of chains. There are analogous hints for a dilute/dense phase 
transition in the case of Stockmayer particles but not for DHS. 

By the end of the 90s there thus was a kind of consensus reached that DHS do not show 
a dilute/dense phase transition because of the competing process of chaining. But then 
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Shelley et al. 



1791 ] perform in 1999 Monte-Carlo simulations for dipolar hard dumbbells 



and spherocylinders. They find a dense/dilute phase coexistence even for almost spherical 
particles and observe, however, increasing sampling problems. They suggest that these 
sampling problems may have prevented the discovery of such a coexistence for DHS. Camp, 
Shelley, and Patey |l80| therefore reinvestigate in 2000 the system of DHS with Monte- 
Carlo simulations. For A = 7.56 they find indeed evidence for o ne o r even two isotropic 
dense/dilute phase transitions. In the same year Camp and Patey |18lj study the evolution 
of a system with 0.001^ < <fi < 0.6^. The simulation shows ring structures for the smallest 4>, 
chain struct ures for intermediate densities (0.06| < < 0.35|), where the phase transitions 
reported in |l8o| | take place, and a s tructurel ess dense liquid for still larger <p. 



Pshenichnikov and Mekhonoshin 



182 



183] also perform Monte-Carlo simulations of DHS 
using a finite system. They find an ordinary dilute/dense transition with a critical point 
already at A c ~ 3 and 4> r = 0. 034. This result is at odds with previous ones possibly 



because the authors of 



1821 . Il83| use instead of the generally preferred Ewald summation 



technique finite cylindrical and spherical systems with 1000 particles. They remark that 



owever, a 



165j | for a 



periodic boundary conditions may cause errors in the case of long-range forces. He 
finite system with N = 2046 particles was previously used by Levesque and Weis 
volume fraction near C a nd m uch higher A without finding a phase separation. 

Morozov and Shliomis jl84| study the intrachain correlations of chains of DHS for large 
A by performing a 1/ A-expansion of the needed statistical integrals. They find a persistence 
length of A/2 for A — > oo. Considering nonideal chains they predict a transition from chain 
to globule structures at A = 3.2 

As discussed above, chaining in monodisperse systems has a strong influence on the 
dilute/dense phase transition. When comparing with experiments done on real ferrofluids it 
should therefore also be noted that polydispersity affects the lengths of chains as theoretically 
predic ted b y Kantorovich and Ivanov 185| and demonstrated in simulations by Wang and 
Holm 12^ . 

The phase behavior of DHS and similar systems i s fur ther complicated when taking into 



account a nonvanishing magnetic field. Sano and Doi 100] investigate the dilute-dense phase 
transition in the presence of a magnetic field in a randomly filled lattice of dipolar particles. 
In zero field they find a rather small critical A and a rather large critical concentration in 
this system. In general a magnetic field makes a phase transition more likely and increases 
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the density differences of coexisting phases. But there is a A-interval near the zero-field A c 
where a phase transiti on occurs for small and large fields but not for intermediate fields. 



Stevens and Grest 



152 



157 



who failed to find a dilute-dense transition for dipolar 
soft spheres in zero field demonstrate the existence of such a transition in Monte-Carlo 
simulations with applied fields. However, the two phases cannot be described as simple 
homogeneous gas and liquid phases since long chains are present in both of them. For both 
dipolar soft spheres and Stockmayer particles they show that A c is getting smaller in applied 
fields whereas the value <p c ~ 0.017 remains mainly unchanged. In |l73| however it is argued 
that this observed phase transition for soft spheres may be just an artifact of the insufficient 
system size. 



Kusalik |l8fij | does not find a phase separation in the system of dipolar soft spheres at 
high fields. However, his value of A ~ 3 is smaller than those considered by Stevens and 
Grest who show that A c > 5 even for infinite fields. The field dependence of A c and <p c in 
he c ase of Stockmayer particles is also studied with Monte-Carlo techniques by Boda et al. 



187]. The aut hors find very similar results to those of Stevens and Grest. Szalai, Chan, 



and Tang [li 



find also qualitatively the same behavior when they study dipolar Yukawa 



particles using the MSA. 

Klapp and Forstmann |189| investigate perfectly aligned DHS, i. e., the infinite-field limit 
within the framework of the RHNC. They also do not find a dilute/dense phase transition 
up to A ~ 5. 

Zubarev and Iskakova investigate an interacting variant of their chain model 



118. 



119, 



and 



that has already been used to calculate the equilibrium magnetization in 
show that in this model chain formation prevents a dilute/dense transition in the infinite 
field limit. 

Morozov and Shliomis 184] find that nonideal chains gain stability against the transition 



to a globular cluster in applied fields. 



B. Ferromagnetic phases 

Here we review papers that have addressed the question whether the classical dipolar 
interaction between the nanoscale particles in a ferrofluid suffices to induce long-range ori- 
entational order of the magnetic dipoles in a liquid state without positional order, i.e., with 
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homogeneously distributed particle positions. 

Already the oldest model for the equilibrium magnetization of dipolar interacting parti- 
cles, namely the Weiss model allows ferromagnetic solutions, i. e., a spontaneously generated 
equilibrium magnetization M^O for H = 0: Eq. (|2.8|) reads in a needle-shaped probe ge- 
ometry for which H = H e 



M = M sat C 



mM 
3k B T 



(6.2) 



when H = H e = 0. This equation has solutions with M^O when M sat mj?)k B T = Xl = 80A 
exceeds 3. Approaching this critical value, xl = 3, of the Langevin susceptibility the initial 
susceptibility (J2.9|) of the Weiss model undergoes a divergence. While the existence of a 
ferromagnetic state for such a rather large range, <p\ > 3/8, of parameters can be ruled 
out experimentally there are experimental hints for a ferromagnetic phase in the paper 
by Mamiya, Nakatani, and Furubayashi from 2000 [19] who find growing ferromagnetic 
fluctuations in their iron-nitride ferrofluid at low temperatures for « 0.15. 

It should be noted that the above mentioned divergence of x disappears when the next 
diagonal term in the (0, A)-expansion, i.e., the 2 A 2 -term is included in the expression for 
the equilibrium magnetization 64]. In this extended Weiss model one has 

"1 20 

— ttpff 

3 9 

where a e g = j^f(H + 4r) is the Weiss mean field. In this case the additional 2 A 2 -term 
prevents the existence of fe rromagnetic solutions for H = 0. 



M = M. 



sat 



' J.2 \2 
-0 A a e ff 



a. 



off 



0(a 5 eS ) 



(6.3) 



Zhang and Widom .192] extend the Weiss mean field theory describing the local field as 
a stochastic variable and find that ferromagnetic phases do not exist for any value of A if 
the volume fraction d> is smaller than some critical value 0.295. 



At the end of the eighties theoretical investigations based on model Hamiltonians 11931 or 

n 

model distribution functions |194j lead to the prediction of ferromagnetic states only for par- 
ticles with sufficient prolateness or oblateness. Nematic ferro- or antiferromagnetic phases 
were investigat ed with different me thods in systems of dipolar ellipsoids and spherocylinders 



for example in 



132 



195 



196 



But in 1992 MD calculations by Wei and Patey |198l Il99| show the existence of fer- 
romagnetic states in the system of dipolar soft spheres for v = 1/1.35 and X/v = 9 
for volume fractions above ~ 0.655. Since both, the soft and the hard sphere poten- 
tial, are not only spherically symmetric but also purely repulsive both should behave very 
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similar. Weis, Levesque, and Zarragoicoechea |200| find indeed in 1992 in Monte-Carlo 
simulations ferromagnetic fluid and solid phases in the system of DHS in their investi gated 
range 0.8? < < 1.2| and A = 6.25. More results for this range are reported 1993 in [137| . 



Weis and Levesque [1621 ] fi nd fe rromagnetic states also at lower volume fractions ~ 0.3? 



Wei, Patey, and Perera 



201] use density-functional methods to investigate the transition 



to ferroelectric phases in systems of dipolar hard and soft spheres reproducing the numerical 



results qualitatively but not quantitatively. More systematic investigations of t 



re phase 



diagrams of Stockmayer particles are undertaken by Groh and Dietrich |202l 12031 |2_04] also 
using density functional methods. The critical temperatures are high. For vq ~ 0.3 for 
example they find a transition between a dilute, isotropic fluid and a dense, ferroelectric 
fluid already for X/vq = 4 o r A ~ 1.3. 

Zhang and Widom [205! investigate the phase diagram starting with the expression for 
the free energy of a van-der-Waals fluid with two additional terms modeling the magnetic 
properties. The authors find qualitatively similar phase diagrams: For high temperatures 
there is a continuous transition from a dilute isotropic to a dense ferromagnetic liquid that 
becomes first order at smaller temperatures. If the short-range interaction energy is small 
enough compared to the dipolar interaction energy then there exists a temperature interval 
where the isotropic/ferromagnetic phase transition is preceded at smaller densities by a 
first-order transition b etwee n two isotropic phases having different densities. 



Levesque and Weis [lfijfl find ferromagnetic states for A = 12.25 and > 0.6?. They 



162 



fo r sma ller was probably 



152 



157] find spontaneously 



however now argue that the ferromagnetic phase reported in 
just a slowly decaying nonequilibrium state. Stevens and Grest 
magnetized fluid phases of Stockmayer particles for A = 4 and > 0.47. For soft spheres they 
find such magnetic fluid phases for the investigated values 4 < A < 9 when is greater than 
a value that depends linearly on A. They see hints for a hysteresis in the magnetic-isotropic 
transition for A = 16. 

Groh and Dietrich |20(t| study the phase transition properties of nonspherical dipolar 
hard particles using a density functional ansatz. In the special case of DHS they find no 
dilute/dense transition between isotropic phases but an isotropic/ferromagnetic transition 
that becomes noncontinuous at A > 1.33 and ~ 0.2. The authors admit that the stability 
of the ferromagnetic phase is overestimated with respect to the simulations. That density 
functional methods are probably not useful for quantitative investigations of the phase be- 
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havior is also poin ted o ut by Osipov, Teixeira, and Telo da Gama 
context by Ivanov ^O^ an d Morozo v 



Klapp and Forstmann 



177. 







207] and in a different 



find fluctuations at the high density boundary of the 



region of existence of isotropic RHNC solutions for the DHS and Stockmayer systems. They 
also find ferromagnetic solutions beyond that boundary. 

Camp and Patey jl81 | report in 2000 that they found in their Monte-Carlo simulations 
ferromagnetic phases at A « 7.5 for <p > 0.57r/6. A Monte-Carlo simulation of a DHS 
fluid with dipolar polydispersity, one with a bidisperse diameter dist ribut ion, and a binary 
mixture of dipolar and neutral hard spheres is undertaken by Cabral |21o| for <fi = 0.8^ and 
A = 6.25. He show s tha t the polydispersity reduces the ferromagnetic order. 

In 2003 Ivanov \lQ^ investigates the BBGKY equation relating the one-particle distribu- 



tion function p(r, m) to the pair correlation function g(l,2). He 



employs 



or the latter the 



self-consistent expression used in density functional approaches [189t I201L I202L |204j . When 
expanding it up to second order in the dipolar interaction the so approximated BBGKY 
equation shows a bifurcation to a solution with spontaneously generated ferrimagnetic long- 
range order for xl > 12(y / 5 — 2) ~ 2.83, i.e., slig htly below the bifurcation threshold for 



ferromagnetism in the Weiss model. Ivanov 



208( | argues that this peculiar ferrimagnetic 



behavior is an artifact that is ultimately due to the mean-field character of the density 
functional ansatz for the pair correlation function. 

Then in 2003 Morozov jsool investi gate s the bifurcation properties of the so-called Lovett- 



Mou-Buff-Gubbins integral equation 



112 



113j | that relates the one-particle probability dis- 



tribution function to the direct correlation function c(l, 2) and that is equivalent to the 
above mentioned first BBGKY equation. He establishes criteria for the bifurcation of solu- 
tions with long-range magnetic order in ellipsoidal and spherical sample shapes in vacuum 
and compares with mean-field predictions, density function theory, and MSA results. Using 
a generalized MSA type expression for c(l, 2) he finds that the strength of short-range cor- 
relations plays a decisive role for the appearance of spontaneous long-range magnetic order 
- the susceptibility diverges when the former exceed a critical strength. More qualitative 
arguments then show that the short-range correlations are most likely to be antiferromag- 
netic. He also shows that approximating c(l,2) by the two lowest orders in its diagram 
expansion always gives rise to liquid ferromagnetic solutions. 

At the end of this section on phase transitions in ferrofluids we should like to quote 
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Teixeira, Tavares and Telo da Gama 



211]. They have presented in 2000 a fairly extensive 



review on the effect of dipolar forces on the structure and thermodynamics of classical 
fluids that covers many aspects of the problems related to gas/liquid condensation phase 
transitions in various dipolar systems. Since their summary of the state of research on these 
problems still gives in our opinion an adequate picture of today's state o f the art — also 
concerning the problem of liquid ferromagnetism — we give two quotes 21 11 ] : "Although 
widely studied in statistical mechanics, the phase diagrams of dipolar fluids in general, and 
of strongly dipolar fluids in particular, have remained largely uncharted" and "In conclusion, 
it is fair to say that a theory is not yet available that is capable of describing dipolar fluids 
over the whole range of densities and dipole strengths. Moreover, the mechanisms driving 
the phase transitions (as well as the location of the phase boundaries) remain unclear." 
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VII. CONCLUSION 

Ferrofluids, i.e., suspensions of magnetic nano-particles in liquid carriers are not only of 
technological interest but they continue to be also the object of many basic research projects. 
Among those has the effect of dipolar interactions that are important in ferrofluids attracted 
considerable interest for two reasons: First, dipolar forces are long-range so that, e.g., the 
equilibrium magnetization (or polarization) is geometry dependent. Second, dipolar forces 
are attractive or repulsive depending on the relative orientation of the particles and their 
dipole moments with the orientational mean being zero. A realistic theoretical description of 
the colloidal suspension will also take into account other interactions between the particles 
and in addition the dispersion of particle diameters and magnetic moments in ferrofluids 
that is absent in molecular systems. The dipolar interaction affects naturally the equilibrium 
magnetization, but also the phase transition behavior. These two subjects were reviewed 
here. 

The question of how dipolar interactions influence the initial susceptibility \ was stud- 
ied for electrically polar systems long before the advent of ferrofluids. The oldest models, 
the Weiss and Onsager models, to capture this influence with mean-field approaches were 
proposed already in the first half of the 20th century. Only after 1970 other theoretical 
methods were used intensively to investigate x by cluster expansions, integral equations, and 
Monte-Carlo and molecular dynamics simulations. Here mainly dipolar hard spheres and 
Stockmayer particles were investigated. While both represent only poorly polar molecules 
like water they are much better approximations to the nanoscale magnetic particles of fer- 
rofluids. In fact the latter seem to be very good realizations of these two model systems. 
Many of the earlier analytical and semi-analytical theories that originally were conceived 
for molecular systems were also extended to account for polydispersity in ferrofluids. 

When comparing cluster expansion methods and integral equation methods with numer- 
ical simulations it seems that the RHNC theory and the simpler perturbation theory by 
Tani et al. j3| agree best with the Monte-Carlo results for Xi & t least for relatively high 
densities. The second theory has also been applied to explain experimental data for a highly 
conc entrated ferrofluid, as has a related cluster expansion result by Ivanov and Kuznetsova 
124|. 

With the synthetization of ferrofluids, i.e., dipolar systems in which saturation could 
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easily be reached not only x but also the full magnetization curve M(H) became of interest. 
Here, the main theoretical results date from 1985 onwards. The proposed theories were 
mostly extension s of those already used to describe \. For example, the one by Ivanov 
and Kuznetsova \yi\ com pares well with experimental data. However, recent numerical 



simulations by Wang et al. 



126, 



128] for x an d M(H) are not well described by this theory 
when the volume fraction of the magnetic particles is small and the dipolar coupling is 
strong: they are better explained by a model of small non-interacting chains of different 
lengths. 

An ever increasing analytical, numerical, and experimental research activity has been 
devoted to phase transition phenomena in dipolar systems, namely, (£) the separation into 
dilute and dense phases and (ii) the appearance of spontaneous (ferromagnetic) long-range 
order in zero magnetic field. Both phenomena refer in the case of ferrofluids solely to the 
subsystem of the magnetic particles without positional long-range order. 

Dilute/dense transitions occur in Lennard- Jones systems and, more generally, whenever 
the attractive part of the interaction is sufficiently strong. In fact simulations have shown 
that switching on an additional dipolar interaction in such systems favors the transition and 
increases the critical temperature. However, the question whether such a phase transition 
occurs also for purely dipolar hard spheres without any isotropic attraction is more difficult 
to answer. The fact that such phase separation has been observed experimentally in real 
ferrofluids does not help much to solve the above theoretical problem when (strong) isotropic 
interactions cannot be excluded to be present in the ferrofluid. 

In the 1970s and 1980s a variety of semi-analytical theories was applied to investigate 
the dilute/dense phase transition behavior of dipolar hard spheres. They all found a usual 
phase diagram with a roughly parabolic coexistence curve. However, there was no agreement 
about the location of the critical point. But Monte-Carlo simulations in the 90s did not 
find such phase transitions for dipolar hard spheres. For strong enough dipolar couplings 
or, equivalently, low enough temperatures a formation of head-to-tail dipolar chains was 
observed instead. They are caused by the highly directional character of the dipolar pair 
potential. Thus, new models were developed to explain the observed behavior by means 
of polymer theory. According to these models the assembly of chains prevent a usual di- 
lute/dense transition or mask it and make its observation in simulations difficult due to 
finite-size effects. 
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Another turning point was reached in 2000 when Camp et al. [180] found evidence for 
one or even two phase transitions betwe en chain- dominated states in Monte-Carlo simula- 



tions. Pshenichnikov and Mekhonoshin 



182L Il83| even observed the common dilute / dense 



transition already at quite low dipolar couplings where chain formation is not very effective. 
They explain the differences to other studies with their choice of a finite system instead of 
periodic boundary conditions. The differences between both simulation methods remain to 
be investigated further. 

The phase transition behavior is complicated further by the presence of a magnetic field. 
In the case of Stockmayer particles the field supports the phase separation, i. e., it enlarges 
the coexistence region. The question whether a field can trigger such a phase separation 
also for dipolar hard spheres is not yet settled. 

In comparison to the question of a dilute/dense transition the research results concerning 
the possibility of having a phase with spontaneous (ferromagnetic) long-range order in zero 
magnetic field seem to be more coherent — albeit only on first sight: Experimental hints for 
such a phase were found by Mamiya et al. 19] and the simulations performed since the 1990s 
generally agree that there exist ferromagnetic phases for dipolar hard spheres and similar 
systems in regions of high densities and strong couplings. This behavior was also found in 
density functional theories although they seem to be less suited to explain the phase behavior 
and the results differ quantitatively. The conditions for the app earan ce of sponta neous long- 
range magnetic order were recently investigated by Morozov |20o| and Ivanov 208 ] using 
integral equation methods. They conclude that mean-field type approximations that are also 
the core of density functional approaches tend to generate artificially long-range magnetic 
order. 

So, all in all we think that a lot of questions/problems related to the equilibrium behavior 
of ferrofluids and dipolarly interacting particles remain to be addressed. 
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FIG. 1: Spheres i and j carrying dipole moments mj and nx,-, respectively, that interact via 
potential ipO)) . 
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FIG. 2: Non-magnetic interaction potential for sterically stabilized ferrofluid particles as function 
of the surface-to-surface distance s. Short dashed line: van der Waals attraction, long dashed line: 
steric repulsion, solid line: combined potential. Note that the scales of both separation and energy 
are different in (a) and (b). 
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FIG. 3: The quantity exp(— v SR ) as function of particle distance for the potential of Fig. |2 (solid 
line). The contribution from the divergent part of the potential at very small surface-to-surface 
distances is suppressed in this plot. The short dashed and long dashed lines show exp(— v SR (r)) 
for best fitting hard sphere and van der Waals potentials, respectively. 
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FIG. 4: Comparison of theoretical results for the initial susceptibility of dipolar hard spheres for 
6 = 0.87T/6: Weiss Q and Onsager Q models, MSA by Wertheim Q, RLHNC results by Patey 
[s^ . RHNC results by Fries and Patey cluster expansion result by Huke and Lucke 6J|, and 
the theory by Tani et al. 26]. Full circles denote Monte-Carlo data taken from a similar figure in 
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